Fermi Acceleration in Plasmoids interacting with Fast Shocks of Reconnection via 

Fractal Reconnection 
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We propose the particle acceleration model coupled with multiple plasmoid ejections in a solar 
flare. Unsteady reconnection produces plasmoids in a current sheet and ejects them out to the fast 
shocks, where particles in a plasmoid are reflected upstream the shock front by magnetic mirror 
effect. As the plasmoid passes through the shock front, the reflection distance becomes shorter and 
shorter driving Fermi acceleration, until it becomes proton Larmor radius. The fractal distribution 
of plasmoids may also have a role in naturally explaining the power-law spectrum in nonthermal 
emissions. 



Introduction - Recent observations of solar X-ray, 
gamma-ray and microwave bursts revealed that energy 
release in a solar flare is very dynamic and that high en- 
ergy particles, i.e. GeV ions and MeV electrons, are gen- 
erated in a very short time period (<1 s). This short time 
variability of bursts, sometimes fractal-like time variabil- 
ity, indicates highly fragmented acceleration regions [l|, 
and these are expected to be above or around the loop- 
top hard X-ray (HXR) source @, Q . HXR spectral ob- 
servations of solar flares have established that efficient 
electron acceleration (10 34 -10 35 electrons/s) occurs dur- 
ing impulsive phase of solar flares Q. To explain the 
high energy particles, several models have been consid- 
ered, such as DC field acceleration inside a current sheet 
[E S 0] , stochastic acceleration in the turbulent recon- 
nection outflow 3 , El, Hoi 1 1 1 and shock acceleration at 
the fast shock [12|, [l3| , though assumed turbulent flows 
are still not revealed. Furthermore, recent studies have 
shown that the role of multiple X-points in a current 
sheet is more important for particle acceleration and en- 
ergy release 14, 15, 

X-ray emitting plasma ejection, shortly plasmoid ejec- 
tion, is one of the direct evidence of magnetic reconnec- 
tion in a solar flare [It}. In the pre- flare phase, a plas- 
moid (a magnetic flux rope) is observed to gradually rise 
up, until it is accelerated upward in the impulsive phase 
in association with HXR burst [3, EH- Multiple plas- 
moid ejection events have been discovered in radio and 
soft X-ray (SXR) observations in association with HXR 
bursts jig . [20( . It is also reported that multiple down- 
flows, which are counterparts of multiple ejections, asso- 
ciate with HXR bursts [21[. These observations may in- 
dicate the relationship between particle acceleration and 
plasmoid ejections. Furthermore, some plasmoids were 
observed falling down to the loop-top in UV and HXR 
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Plasmoid ejection drives dynamic feature of magnetic 
reconnection. It is known that plasmoid ejection in- 
duces inflow into the current sheet and increases recon- 
nection rate 24, [5^]. Plasmoid formation repeats self- 
similarly in a current sheet and makes fractal-like or tur- 



bulent structure in a current sheet via fractal reconnec- 
tion [H, [H, E3, El HI El EH, S3 • Plasmoids in a current 
sheet are unstable for coalescence instability and repeat 
lots of collisions with each other, then merging to a single 
large plasmoid which is finally ejected out of the current 
sheet. During this process, strong DC-field is enhanced 
between the two colliding plasmoids and shrinkage of a 
large plasmoid also accelerates particles impulsively, un- 
til a lar g e st p lasmoid is ejected outside of the current 



sheet [141 1331 . |34j, [35J, |36(. The large numbers of plas- 



moids remain as the exhaust impacts the loop-top shock. 
This is directly observed as multiple plasmoid e ject ions 
and downflows correlated with HXR emission [2fj, |21 | 
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FIG. 1. Particle acceleration driven by multiple plasmoid 
ejections colliding with the fast shock. Multiple plasmoids of 
various scales are intermittently ejected upward and down- 
ward out of a turbulent current sheet and collide with the 
termination shocks of reconnection outflows, i.e. fast shocks 
above and below a reconnection X-point, where particles are 
effectively accelerated via shock acceleration process trapped 
in a plasmoid. 
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and indirectly as power-law distributions of HXR and 
U V footpoint bright points 0, |37} . 
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FIG. 2. Scenario of shock acceleration at the fast shock 
trapped in a plasmoid; particles trapped in a plasmoid col- 
lide with the fast shock, when they repeat reflections up- 
stream of the shock because of magnetic pressure gradient 
(mirror force). During the passage of a plasmoid through the 
shock front, trapping distance of particles becomes shorter 
and shorter and drives Fermi acceleration process, until it be- 
comes microscopic scales enough for particles to escape from 
trapping. 

In this paper, we propose acceleration model 
during the dynamic process of interaction between 
bi-directionally ejected fractal plasmoids and fast shocks 
just below and above the current sheet. Here we focus 
on the fractal distribution of plasmoids and their role 
for the trapping favorable for acceleration. 

Scenario of First-order Fermi Acceleration - We pro- 
pose a Fermi acceleration process when multiple plas- 
moids collide with a fast shock above the loop-top. Once 
reconnection occurs in the corona, reconnection outflow 
generates termination shocks (fast shocks) above the 
loop-top and below the magnetic flux rope [38j ]. The 
scenario we propose is as follows: (1) multiple plasmoids 
formed in a fractal current sheet are ejected downward 
(upward) and collide with a fast shock below (above) the 
current sheet with the speed less than or comparable to 
the Alfven velocity (300-1000 km/s) with trapped parti- 
cles inside (Fig. 1). (2) During the collision, magnetic 
pressure gradient at the shock front, i.e. magnetic mirror 
force, reflects and traps particles upstream the fast shock 
(Fig. 2). (3) Through the passage of the plasmoid, the 
distance between two reflection points along the field line, 
L, becomes shorter and shorter. At that time, each par- 
ticle gets momentum from the shock front A(mv)=2mu 
via one collision, resulting into the first-order Fermi ac- 
celeration. Here m and u are particle mass and the rel- 
ative velocity of the shock front to the rest frame of the 
upstream. (4) Finally when the bouncing timescale be- 
comes comparable to the reflection time scale determined 
by ion cyclotron frequency (2-10)S7~ 1 [H, |4(|, in other 
words, when the distance L becomes comparable to the 
Larmor radius, particles escape from the trapping up- 
stream the fast shock. 

During this process, plasmoids disturb the ambient 



plasma and scatter particles. It is also known that the 
parallel motion of a particle tends to change into the per- 
pendicular one at the rotational discontinuity. The same 
process may happen when magnetic field lines are bent 
at the fast shocks. Since the pitch angle increases dur- 
ing the passage of the shock front, it would positively 
work for the trapping of particles. Similarly, betatron 
acceleration by compression of a plasmoid downstream 
and even upstream the shock will play the same role in 
continuing reflection longer. The guide field enables par- 
ticles to move parallel to the shock and further DC-field 
acceleration occurs. 

For the first-order Fermi process, e nerg y gain of parti- 
cles is described by dE/dt = 4uE/L [4l|. Particles con- 
serve the action parallel to the field lines, so that the time 
variation of particle energy is E(t) = C/L 2 =C/(Ln-2ut) 2 , 
where C is constant. This leads to the ratio of acceler- 
ated and injected particle energies written by the ratio 
of Larmor radius of accelerated particles i'L,acc and the 
size of a plasmoid Lq, 




Here the adiabatic invariant leads to E Lq= E acc r 2 L acc = 
E acc { 1 mcleB) 2 =E acc {E acc /eB C ) 2 =ElJ{eBc) 2 . Then, 
we finally get the following equation, 

- 4640 (^Y f-^-V* (2) 

E Q 4b4U Vl0 4 nJ [lOGj \WkeVj () 

Here we note that Larmor radius in the corona is 1- 
100 m, so that the scale gap between Larmor radius and 
the size of a plasmoid (10 4 -10 8 m) quite positively works 
on this acceleration process. If particles fall into the loss 
cone before the end of the bounce motion, the energy 
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FIG. 3. Particle energy spectrum calculated by solving equa- 
tions of motion of test particles, where dotted and solid lines 
show energy spectra before and after acceleration, respec- 
tively. We assumed a cylindrical flux tube and compressed 
it below the shock front with the compression ratio less than 
4 moving upward with Alfven velocity in the rest frame of the 
upstream, in which particles are trapped and accelerated by 
the electric field inside the shock front of finite width. 
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gain would be determined by the final loop length. The 
time scale of acceleration corresponds to the transit time 
of a plasmoid through the shock front, 
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Tacc = = 10 
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(3) 



This means effective acceleration in a short time period. 

This acceleration process is without the selection of 
protons and electrons. Electrons gain relatively small 
energy compared with protons but complement by the 
number of collisions. The injection energy necessary for 
this acceleration process is determined by the two con- 
straints. One is that initial particle velocity is supcr- 
Alfvenic, when a trapped particle can repeat bouncing 
upstream the fast shock. The other one is that acceler- 
ated particles overcome the energy loss rate by Coulomb 
collision fcgT/T e i, where ks is Boltzmann constant and 
r e i is collision time between electrons and protons such 

that f = f-4.9x ICr 9 (^) (j^)-* [keV/s] > 
for electrons, where E is electron energy and n is elec- 
tron density. This constraint leads to the lower cut off 

energy of electrons, E c = 13 ( 10lu ^ m _ 3 ) 1 (f) 1 [keV] 0. 

For protons, E c = 0.11 ( 10 ii£ m - M )* [MeV]. These 

constraints are greater than super- Alfvenic and require 
additional acceleration in the initial phase. Electrons 
are heated up by slow shocks elongating from the recon- 
nection point, from Eq to Eq//3 where f3 is plasma beta 
(^0.01 in the corona). On the other hand, protons could 
be accelerated in multiple X-points via merging of plas- 



moids 13J,l35|,|36| i.e. E = eEd=ev in Bd= 100 keV [v in /l 
km s _1 )(.B/100 G)(g?/10 km) as well as Fermi process in 
a shrinking plasmoid after the coalescence [13, HH , where 
d is acceleration distance and Vi n is inflow velocity. 

It is interesting to note that multiple internal shocks 
can be generated in reconnection outflows during the 
nonlinear evolution with finite fluctuation at the dif- 
fusion region [42[. At that time, multiple small-scale 
plasmoids can interact with the internal shocks in a 
fractal manner as well as large plasmoids do with the 
loop-top shocks. 

Integrated Power-law Spectrum — The energy spectrum 
of particles accelerated in a plasmoid colliding with a fast 
shock is derived from the equation of mass conservation 
in energy space, 

lH + 0E{-0l N ^ t) )=-- (4) 

where r is the run-away time scale much larger than the 
acceleration time t acc and hereafter we neglect the first 
term dN/dt (if we consider time dependence of N(E,t), 
a factor exp(-at/V) is multiplied to the steady solution). 
For the first-order Fermi process, dE/dt = AuE/L = 
iuE 3 / 2 /C 1 / 2 because EL 2 =C. If we assume the relative 



shock velocity along the field line u is independent of 
particle energy E, we get the solution 



N(E) = E- 3/2 



exp 



G 



2utVE 



oc 



(5) 



However, the relative shock velocity u actually de- 
pends on the distance of bounce motion along the closed 
field lines between the two reflection points of a plas- 
moid L(E), such that u(E) = Uo/ sin | = uo/ sia(L/2R)— 
uq/ sin(v / C /2R-\/E), where uq is the relative velocity of 
the shock front normal to the rest frame of upstream 
(~Alfven velocity) and 8 is the central angle of the cur- 
vature L (9=L/R; R is the radius of curvature). Then 
we get the solution, 



N(E) = £T 3/2 exp 



C 



2utE ! 2 



du/dE 



dE . (6) 



In an extreme case (0 <Sl), the relative shock velocity is 



log N(E) 



N(E) oc E"P' 2 



slope: -p/2 




FIG. 4. Superposed power-law distribution of multi energy 
spectra of plasmoids colliding with the fast shock. 

u(E) = 2u /9 =2u Q R/L =2^VE, therefore, du/dE = 
u/2E. This simplifies the equation (6) to 



N{E) 



E 2 exp 



C 



8u tRE 



oc E~ 



(7) 



which asymptotically approaches to E~ 2 for E ^>1. 

For comparison, we performed test particle simulation. 
We assumed 2D magnetic field configuration of a closed 
magnetic island colliding with a fast shock with ana- 
lytical model whose downstream side is constantly com- 
pressed in the rest frame of a plasmoid, in which trajec- 
tories and energies of protons were calculated by solving 
gyromotion. The initial energy spectrum is soft power- 
law distribution with power-law index 7 (Fig. 3). This 
numerical simulation shows that particles are effectively 
accelerated at the shock front and the power-law spec- 
trum grows harder to the power-law index -2. This is 
consistent with the previous estimation (detailed discus- 
sion in Nishizuka ct al. in prep). 
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The observed HXR spectrum would be the superposi- 
tion of several HXR spectra emitted from the numbers 
of plasmoids colliding with a fast shock. Both acceler- 
ated particle energy and acceleration time are propor- 
tional to the size of a plasmoid Lo, as shown in equa- 
tions (l)-(3). Therefore, if the size of a plasmoid Lo 
depends on the power-law distribution via fractal recon- 
nection process, the power-law distribution of the ob- 
served HXR spectrum can be naturally explained, even 
if some energy spectra from plasmoids does not develop 
to the power-law distribution enough. Each spectrum 
tends to have the peak at the temperature determined by 
the adiabatic heating, i.e. E = (Lq/tl^Eq. Then, here 
we assume the power-law distribution of plasmoids with 
power-law index a, i.e. N(Lq) cx L^ a (fractal current 
sheet), and the power- law energy spectrum with power- 
law index p, N(E) oc E~ p . At that time, these two 
power-law indices a and p are related by the equation 
N(L )dL = N(E)dE, simplified to p = \{a + 1). With 
numerical simulation, the power-law index a is shown to 
be a=2 [3l], [37J • If we adopt a = 2 based on the simula- 
tion result, this equation gives us p — |, which is harder 
than the observation, though it can become softer if a 
varies. This may also indicate the possibility that we can 
expect microscopic parameter a from the observed HXR 
spectrum. 

Our model is a unified model of magnetic reconnection 
and particle acceleration in a solar flare, in which 
particles are effectively accelerated coupled with the dy- 
namics of plasmoid ejections colliding with fast shocks. 
Impulsive ejections of fractal plasmoids explain observed 
short time variations of HXR and 7-ray bursts (<1 
s), HXR power- law spectrum and maximum energy of 
acceleration (MeV electrons and GeV protons). Multiple 
ejections of plasmoids also generate turbulent outflows 
and internal multiple shocks (42| . These would further 
enable stochastic acceleration. The trapping between 
fast shocks and plasmoids increases local coronal density 
and may explain observed HXR flux via bremsstrahlung. 
Our model can be applied to fast shocks both at the 
loop-top and below the flux rope, which is consistent 
with double coronal HXR sources, and then also applied 
to astrophysical jets with high energy particles. 

We acknowledge the supports in part by a JSPS Re- 
search Grant, in part by the JSPS Core-to-Core Program 
22001 and by the Grant-in-Aid for the Global COE Pro- 
gram "The Next Generation of Physics, Spun from Uni- 
versality and Emergence" . The numerical simulation was 
carried out on Altix2700 at YITP in Kyoto University. 
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